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Q , Abstract 

\o . 

Q\ • We present a lattice investigation, partially including the effects of dy- 

namical quarks, of the heavy-light mesons using NRQCD for the heavy 
quark and the Wilson action for the light quark. We performed an exten- 
p) • sive calculation of the spectrum employing a multi-state, multi-exponential 

O . fitting analysis which enabled us to extract the 2S—1S splitting as well as 

the ^Pi—S and hyperfine splittings. The Wilson action introduces a large 
systematic error into the calculation, and within this uncertainty we obtain 
^ ■ agreement with experiment. We performed a comprehensive calculation of 

the heavy-light meson mass, investigating three methods and their range 
of validity. The agreement found between these methods confirms that 
Lorentz invariance can be restored at this order in NRQCD by a constant 
shift to all masses. We calculated spectroscopic quantities over a wide 
range in heavy quark mass and were able to perform a detailed investi- 
gation of heavy quark symmetry around the b quark mass. In particular, 
we extracted the nonperturbative coefficients of terms in the heavy quark 
expansion of the meson binding energy. We demonstrate the importance 
of using tadpole-improved operators in such a calculation. 
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1 Introduction 



Lattice field theory has an important role to play in the study of 5-Physics. 
Through this approach it is possible not only to perform a first principles cal- 
culation of experimentally important quantities, but also to test heavy quark 
symmetry predictions of Heavy Quark Effective Theory (HQET) and extract the 
nonperturbative coefficients of the heavy quark expansions which are required to 
extract numerical results from HQET. Furthermore, comparison with the pre- 
dictions of various models for these coefficients provides a test of the underlying 
assumptions of the models. 

The first stage of a lattice study of B mesons is the calculation of the spectrum. 
As the simplest quantity to extract, it is the foundation for more complicated 
calculations of B decay processes. A correct reproduction of the spectrum gives 
confidence in the calculational procedure, as well as providing predictions for 
states not yet measured experimentally. In the spectrum, several physical scales 
are important. The gross structure of the B spectrum is due to excitations of the 
light quark and hence mass splittings related to these, such as P—S and 2S—1S 
are of 0{Aqcd)- Conversely, the heavy quark is approximately a static colour 
source and its spin affects the spectrum at the scale of the hyperfine splitting, 
0{AlaD/M). 

Thus, to reproduce the spectrum, the systematic errors arising from both the 
heavy and light quarks must be rigorously controlled. Previously, simulating the 
b quark (as opposed to a static quark) on the lattice presented a major problem 
since the 0{Ma) discretisation errors of the standard relativistic fermionic action 
were 0(1) for a b quark on a typical lattice, ~ 2.0 GeV. However, the 
success of the non-relativistic QCD (NRQCD) approach in calculations of T 
spectroscopy ^ shows that the b quark can be simulated directly on lattices 
currently available. In NRQCD the heavy quark mass scale is removed and 
the action becomes an expansion in the typical (non-relativistic) heavy quark 
velocity, which is a small parameter. The precision of the calculation is controlled 
simply by the order at which the expansion is truncated. For heavy-light mesons, 
the terms included in the action at each order in the heavy quark velocity is 
determined by the order of the terms in 1/M; this is discussed in more detail 
in reference [^. Considering the almost spectator role of the heavy quark in the 
heavy-light meson, we include only first order i.e. 0{1/M) terms in the NRQCD 
action in this study. 

For light quarks, recent studies have shown that the 0{ma) discretisation 
errors of the Wilson action are significant for a high statistics calculation on 
typical lattices even for small bare masses (i.e. the constituent mass seems to 
be the important scale for 0{ma) effects). Improved light fermion actions, for 
example the tadpole- improved Clover action with errors of roughly O(m^a^), are 
needed. However, in this initial calculation we have employed the Wilson action 
and expect this to give rise to the dominant systematic error. 
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Another uncertainty in connecting lattice calculations to the physical world 
is the quenched approximation. In this study we have partially included the 
effects of dynamical quarks. In addition, these computations are performed in 
conjunction with a simulation using (3 = 6.0 quenched configurations 0] and in 
the future we expect to be able to isolate and remove the full effects of quenching. 

The organisation of this paper is as follows. The details of the simulation 
are given in section ^ followed by a discussion of the expected systematic errors. 
Particular attention has been paid in our analysis to the fitting procedure and 
extraction of ground and excited state energies. Our methods are outlined in 
section | and illustrated with a subset of the data. We use the results over a range 
of heavy quark masses to investigate heavy quark symmetry in the meson binding 
energy in section p.2| . The lowest order coefficients of the heavy quark expansion 
A and < —D^ > are extracted. The heavy-light meson mass is calculated in 
section |5.2.1| ; three methods are compared and their range of validity investigated. 
In section we turn to the mass splittings, B*—B, B**—B and Bg—B, and 
compare the simulation results with experiment and the expectations of heavy 
quark symmetry. The conclusions and a summary of the spectrum are given in 
section |. 



2 Simulation Details 

The simulations were performed using 100 16'^x32 gauge configurations a.t (3 = 5.6 
with two flavours of staggered dynamical sea quarks with a bare quark mass 
of amsea = 0.01. These configurations were generously made available by the 
HEMCGC collaboration; more details can be found in |Q. We fixed the con- 
figurations to the Coulomb gauge. The light quark propagators were generated 
using the Wilson fermion action, without an 0(a) improvement term, at two val- 
ues of the hopping parameter, k = 0.1585 and 0.1600. The former corresponds 
to a quark mass close to strange, where Kg = 0.1577 from M^, while 0.1600 is 
somewhat lighter, with = 0.1610. 

In this simulation we truncate the NRQCD series at O(l/Mo) and the action 
takes the form: 

S = Q^{Dt + Ho + 6H)Q (1) 

where 

A(2) a ■ B 

Ho = — and 6H = -cb^^. (2) 

2Mo 2Mo ^ ^ 

Tadpole improvement of the gauge links is used throughout, where Uq = 0.876 
measured from the plaquette, and the hyperfine coefficient is given the tree-level 
value cb = ^- We use the standard Clover- leaf operator for the B field. A constant 
mass term has been omitted from the action for simplicity of the calculation. It 
has no effect on mass splittings and amplitudes and a suitable energy shift for all 
masses can be added post-simulation. This will be discussed in further detail in 
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section [5.2. 1| . The heavy quark propagators were computed using the evolution 
equation: 



In the static hmit this reduces to 

Gt+1 - UlGt = 5.,o (4) 

The parameter n is introduced to stabihse unphysical higher momentum modes 
as Mo is decreased and from perturbation theory one estimates n^3/Mo. This 
is only a rough guide and for a high statistics simulation a more conservative 
lower limit on n may be needed to ensure a smooth transition between regions of 



different n as Mq decreases. This point will be illustrated in Section [O . 

The heavy quark propagators were computed over a range of values of Mq to 
probe both the lower limit of NRQCD where the perturbative couplings diverge 
and the theory breaks down, around aMo^O.6 and the aMo — > oo limit where 
signal / noise becomes poor as NRQCD tends to the static limit. Note that for 
heavy-light mesons there is no upper limit imposed on Mq by discretisation errors 
in the heavy quark motion, apq ~ aAgcn ~ 0.2 <^ 1 independent of Mq (where 
we take aKqcn = oAy = 0.185 for these configurations). This is in contrast to 
heavy-heavy mesons where the kinetic energy, not the momentum, of the heavy 



quark is approximately Aqcd and apq ~ ^o^MAqcd- Thus, discretisation errors 
naively become large when apq ~ 1 i.e. aM ~ l/aAqcn ~ 5. We generated 
heavy quark propagators at 11 values of (aMo,n) corresponding to (0.8,5), (1.0,4), 
(1.2,3), (1.7,2), (2.0,2), (2.5,2), (3.0,2), (3.5,1), (4.0,1), (7.0,1) and (10.0,1), and 
the static limit. This roughly corresponds to a range of meson masses from Ms/2 
to 4Mb and is sufficient for a reasonable investigation of heavy quark symmetry. 
However, it is not possible to simulate the D meson on this lattice using NRQCD; 
a larger lattice spacing is required, or /3"-''=°^5.85. 

A NRQCD heavy quark propagator describes the forward propagation of a 
quark. To increase statistics we applied the time reversal transformation to both 
the gauge fields and the light quark propagators and repeated the heavy quark 
evolution; this is equivalent to calculating the propagation of the backwardly 
moving anti-particle. 

Following the methods of |1|, |^ meson correlation functions corresponding to 
both S and P states were constructed from the quark propagators using interpo- 
lating operators at the source of the form: 

o=Y. Q^rt(xi-f2)g(x2), (5) 

Xl,X2 
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n 


'^0 (0-+) 


I 


(1-) 




'Pi 


A, 


^Po (0++) 




(1++) 


AjCTj — AjfTj 


(2++) 


Aj(Tj — AjfTj 




AjCTj + AjCi 







Table 1: The lattice operators, fi, and associated quantum numbers used in the 
simulation. f2 is a 4 x 2 matrix in spin space. Note that since charge conjugation 
is not a good quantum number for the heavy-light system there is mixing between 
the ^Pi and ^Pi states. 



and at the sink, where F = Q(j){xi —X2) specifies the quantum numbers {Q) of 
the state, detailed in table |l|, and the smearing (0) applied to the heavy quark. 
Smearing functions motivated by the hydrogen model were chosen to project out 
the ground and first excited state; for S-states 

01 (r) = e-"^/"" 

02(r) = (1 - r/(2ro))e-'^/(2ro) (g) 

where tq = 3.0. Similarly for P-states 

01 (r) = (-1 -r/(2ro))e-''/('"») 

02(r) = (r/(3ro))'e-^/(3'-") (7) 

It is sufficient to use the same radius, ro, for all Mq since the heavy quark is 
almost a spectator in the heavy-light system and the meson wavefunction is not 
expected to change significantly with Mq. A delta function was also used for both 
S and P states. 

With the aim of performing a multi-state, multi-correlation function analy- 
sis the heavy quark was smeared using all the source-sink combinations of the 
smearing functions. The meson propagator is then 

yi,y2 

where 

Gtivi) = T.Gtiy- x)r(^'^)(f)e^i^ (9) 
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and sc, sk = 1,2,/ corresponding to the ground state, first excited state and 
delta function smearing functions respectively. is the light quark propaga- 
tor. Finite momentum ^5*0 correlators were generated for a subset of the data, 
aMo = 0.8—4.0, at a\p\ = 1 and \/2 in units of j^. To improve statistics, where 
appropriate, the correlators were averaged over different polarisations and differ- 
ent momentum directions. 

3 Systematic Errors 

There are systematic errors associated with both the heavy and the light quark 
in the heavy-light meson. The truncation of the NRQCD series at 0(1/ M) in- 
troduces an absolute systematic error of 0{Aqcd{-^qcd/MY) associated with 
the heavy quark. This corresponds to a relative error of 0{Aqcd/M) in quan- 
tities of 0{Aqcd{Aqcd/M)), naively estimated to be 0{v) ~ 10%, and an 
0((Aq(7£)/M)^) ~ 1% error in quantities of 0{Aqcd)- The use of tree-level 
coefficients for the relativistic corrections will introduce errors of the form 
0{asAQCD{AQCD/ M)) , which are of the same magnitude as those above for as ~ 
0.1. Similarly, the magnitude of the scaling violations due to using the Wilson 
action for the light quark are O^Aq^^a), or a relative error of O^Aqcdo.) ~ 20% 
for this ensemble. Hence, the systematic errors associated with the light quark 
are expected to dominate and the finite volume errors for a lattice of size ~ 1.6 fm 
are unlikely to be larger than this for a B meson. It is important to note that this 
is a high statistics calculation with statistical errors of around 5%, much smaller 
than the systematic errors. 

An indication of the size of the light quark uncertainties in the simulation 
can be found by comparing the lattice spacing obtained from various quantities. 
Table ^ details the values of extracted from various light spectroscopic quan- 
tities for this ensemble, where the Wilson action was used for the light quarks. 
Also included is the inverse lattice spacing obtained from the T spectrum using 
NRQCD |jl|. In this case all terms in the NRQCD series were included up to 
0{Mv^) and the estimated systematic error is ~ 1%. Clearly there is marked 
lack of agreement in the values obtained from observables dominated by presum- 
ably similar typical physical scales indicating large light quark uncertainties. In 
addition there is disagreement between the scale obtained from light and heavy 
spectroscopy. With partial inclusion of quark loops at finite sea quark mass, the 
running of the strong coupling does not match that in the real world but should 
provide an improvement on the quenched approximation. At present the light 
quark uncertainties in the spectroscopy hide any indication of this and the vari- 
ation in the lattice spacing is similar to that found for a lattice with a similar 
(/? ~ 6.0) in the quenched approximation. Currently, the light quark spec- 
troscopy is being re-analysed using a multi-smearing, multi-state analysis with 
Clover fermions [Q. 
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A large uncertainty in translates into similar uncertainties in predictions 
for physical quantities. While it is useful to quantify these by using a range of 
values for to convert into physical units, observables will exhibit systematic 
errors differently and the resulting uncertainty may not necessarily be expressible 
by merely using a range of values for a~^. With this caveat we will take — 
1.8 — 2.4 GeV. The emphasis of our analysis will be on determining the heavy 
quark mass dependence of various quantities and the nonperturbative coefficients 
of the 0(1/M) terms. 

With light quark propagators at only two values of ki and with a much worse 
signal to noise ratio for the data at ki — 0.1600 compared to ki — 0.1585, it is 
not possible to perform a trustworthy chiral extrapolation. In addition, the light 
quark uncertainties affect the determination of the bare strange quark mass. Kg 
extracted from pseudoscalar mesons, using the lowest order chiral mass depen- 
dence to find the 'experimental mass' of the pure ss pseudoscalar, differs from Ks 
obtained from the vector meson 0. Thus, most of the results presented will be 
for Ki = 0.1585. 



"force" 


2055 




2140 


Tfip 


1800 


T(NRQCD) 


2400 



Table 2: a ^ in MeV obtained from various observables for this ensemble. The 
"force" is calculated from the static quark potential using tq. Statistical errors 
are 50 - 100 MeV. 



4 Analysis 

With three smearing functions for the heavy quark in all source-sink combina- 
tions we generated a 3 x 3 matrix of meson propagators, C.ij(t), for each state, 
where i, j = 1, 2, /. With several independent determinations of the meson corre- 
lation function it is possible to perform simultaneous multi-exponential fits and 
thus rigorously constrain the ground state energy and amplitudes of the meson; 
higher excited states can also be obtained. Using a subset of the data involving 
correlation functions smeared at the source and local at the sink, vector fits were 
performed using the fitting function, 

Cu{t) = Y.h{i,n)e-^-' i = l,2 (10) 

n=l 
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where Ugxp is the number of exponentials used in the fit. Similarly, 

Tlexp 

Cij{t) = E a(i,n)a(j»e-^"* ij = 1,2. (11) 

n=l 

was used to fit to the 2x2 matrix of correlators formed from the correlators 
smeared using the ground and first excited state smearing function. 

The fitting procedure applied to extract the mass spectrum from these corre- 
lation functions can be expressed as follows: 

1. Choose the end of the fitting range, tmax, to be a timeslice where the 
noise/ signal^S. 

2. Begin with a 1 exponential correlated fit. 

3. Vary the beginning of the time range, tmin, over all possible values i.e. 
ncoritmax — tmin + 1) ~ ^par > 0, where Ticor is the uumber of correlation 
functions and Upar the number of parameters in the fit. 

4. Accept fits for which the quality of fit parameter Q > 0.1. 

5. In the case of multi-correlation function fits repeat (3) including 2 then 3 
exponentials. 

For each set of fits with a fixed number of exponentials, the criteria for the 'best' 
fit is taken to be the earliest value of tmin where Q has reached a plateau i.e. is 
stable to successive increases in tmin- 

With each addition of another exponential to the fit function the range of 
tmin with an acceptable Q moves closer to the origin and the extra exponential 
removes the excited state contamination from the lower state. A plateau is sought 
in the ground state fit parameters obtained using 1, 2 and 3 exponential fits. 
To achieve this and resolve the first, and higher, excited states there must be 
a sufficiently different overlap with the ground and excited states between the 
correlation functions with different smearing combinations. In addition, a signal 
for isolating the first excited state must be a non-zero overlap of the state with 
the excited state smearing function. In general, the highest excited state is likely 
to be untrustworthy, i.e. a 3 exponential fit is needed to provide confidence in 
the first excited state. 

In order to perform correlated fits to quantities as a function of the meson 
mass 100 bootstrap ensembles were generated for each correlated fit to the prop- 
agators. However, we found that our statistics were insufficient to be able to 
perform the bootstrap procedure above Uexp = 1- Thus, only Uexp = 1 fits were 
used. It was possible to do this with confidence only because we were able to 
check agreement with the multi-state, multi-smearing analysis of the unbooted 
correlation functions. 
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In the construction of the correlation functions we are combining hght quark 
propagators which include the contribution of the backwardly propagating anti- 
quark with the NRQCD heavy quark propagating solely in the forward direction. 
This leads to a unphysical contribution to the meson correlation function which 
should be small at least through tmax ^ Lt/2 = 16. Examination of the correla- 
tion functions indicates the contribution is small for tmax^'24:. 

5 Results 

5.1 Fitting results for Esim 

The simplest quantities to extract are the ^5*0 and ^Si energies, E^^^ and -E]^^ 
respectively, which determine the exponential fall-off of the S-state correlators. 
With the omission of the constant mass term in the NRQCD action Egim is not 
the meson mass, M2; the calculation of this quantity is discussed in the next 
section. However, physical mass splittings can be obtained from differences in 
Esim for different mesons. In particular, to have confidence in the results it is 
vital to ensure the excited state contamination of Esim is minimal. We use the 
aMo = 1.0 and ki = 0.1585 results to illustrate how we achieve this with our 
fitting analysis. 

Figure shows the effective masses of all the ^5*0 smeared correlation func- 
tions. The correlators smeared with the excited state smearing function have a 
larger contribution from excited states compared to those smeared only with the 
ground state smearing function, as required to separate ground state and excited 
states. Increasing Mq from 1.0 to 10.0 corresponds to an approximate 7-fold in- 
crease in the meson mass. However, a comparison with the effective mass of Cu 
for aMo = 10.0 in figure |^ shows there is only a gradual change in the meson 
wavefunction, as expected from a hydrogen-like picture where the heavy quark 
plays a minor role. 

The ^Sq ground state energy extracted from the Uexp =1,2 and 3 vector fits 
is presented in figure |^ as a function of tmin, with tmax fixed at 20. Fits are shown 
for which Q > 0.1 and for which all states included in the fit are resolved. For 
example the riexp = 2 vector fit provides results for Esim consistent with n^xp = 1 
for tmin > 10. The first excited state, however, is no longer resolved as seen by 
an excited state amplitude consistent with zero. A ne^p = 3 fit can be performed 
using tmin = 2 and 1. This does not give sufficient additional timeslices, however, 
to resolve the third excited state. 

A noticeable feature is the 2cr 'wiggle' in Esim which also appears in the 
effective mass of the Cuit) correlator. There are two possibilities, excited state 
contamination or finite statistics. If the former is responsible the contribution 
from excited states is small enough that it cannot be resolved using Uexp = 2. 
Thus, we quote the uncertainty in Egim due to the fitting procedure to be 2a. 
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Figure 1: The effective masses of the ^5*0 meson correlators for aMo — 1.0 and 
Ki = 0.1585. 
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Data 




fit range 


Q 


n E ■ 


oE' ■ 

Sim 


11,21 


2 


3-20 


0.2 


0.484(2) 


0.83(2) 


11,12,21,22 


2 


3-16 


0.8 


0.485(2) 


0.80(2) 


11 


1 


3-20 


0.9 


0.485(3) 





Table 3: A comparison of the energies extracted from fits to 3 data sets of So 
correlators for aMp = 1.0 and ki = 0.1585 



This feature is seen in the results for all values of Mq. 

We found -E^j^ ^^e first excited state extracted using riexp = 2 to be 
stable around 0.8 and the corresponding amplitude to be well determined for 
tmin = 2—6; in this region there is a significant excited state contribution to C21 
while the ground state dominates Cu. A Uexp = 3 fit over a wider range in tmin 
is needed to provide confidence in this result. Cu was not found to be useful in 
providing a larger overlap with the first or second excited state since even higher 
excited states were also present with significant contributions. 

A similar picture is found in the results for the matrix fits presented in figure ^ 
Esim and E'^-^ are consistent as tmin and n^xp are varied and in agreement with 
those obtained using the vector fit. This is proof that we have isolated the ground 
state and minimised excited state contributions to Egim- 

As mentioned previously, we were limited to bootstrapping only ng^p = 1 fits. 
So we also performed Ue^p = 1 correlated fits to the Cii(t) correlation functions. 
The effective mass and values of Egim obtained, with the corresponding values of 
Q, are shown in figure ^. A plateau can be seen in the variation of Egim with tmin 
from tmin = 3. Table ^ compares the best fits from all three fitting methods. The 
values for the ground state using multi-exponential fits are consistent with the 
nexp = 1 fit to Cii(t). The statistical error in the fits to Cn is larger. However, 
this will not affect quantitative results of the analysis as the systematic errors 
are much larger than the statistical errors. 

The fitting range of 3—20 was found to be optimal for the ^5*0 state for all Mq 
at Ki = 0.1585, while for ^5*1, 4—20 was used. Table H details the corresponding 
values of E^^^ and E^im- The first excited state energy extracted using Uexp = 2 
correlated fits to the vector of smearing functions with the fitting range 3—20 for 
qMq = 1.0, 2.0 and 4.0 are presented in table |^. The E'^-^—Egim splitting is of 
the order of qAqcd ~ 0.2, and only weakly dependent on the heavy quark mass, 
as expected for the excitation of a light quark within the heavy-light meson. 

We were not able to perform such an extensive analysis in the static limit. 
Figure 1^ presents the effective masses for the Cu, C21 and Cu correlation functions 
for a static heavy quark and = 0.1585. The static results suffer from the signal 
disappearing around timeslice 11. The smeared correlators C12, C21 and C22 are 
dominated by noise and are not useful. We performed rtexp = 1 correlated fits to 
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Figure 2: The variation of ^5*0 ground and first excited state energies with tmin 
extracted from vector fits to the Cu and C21 meson correlators for oMq = 1.0 
and Ki — 0.1585. Fits performed with Uexp = 3 are shown as squares, diamonds 
represent Uexp = 2 and crosses Uexp = 1- Only fits for which Q > 0.1 are shown. 
The variation of Q with tmin is shown as the solid line for Hexp = 1 fits, the 
dotted line for Uexp = 2 and the dashed line for rigxp = 3. tmax is fixed at 20. Also 
presented is 6(2, 2), the overlap of the first excited state with the C21 correlator. 
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Figure 3: The variation of ^5*0 ground and first excited state energies with i^j„ 
extracted from matrix fits to the Cn, C12, C21 and C22 meson correlators for 

aMf) = 1.0 and k,| = 0.1585, where t^ax = 16. The same symbols are used as in 
the previous figure. Also presented is the overlap a(2, 2) of the first excited state 
with the excited state smearing function. 



12 



0.47 



5 10 
t„,„ 



0.8 



0.4 




0.0 



5 10 
t„,„ 



Figure 4: The variation of E^^^ with tmin extracted from a single exponential fit 
to Cii for gMq = 1.0 and ki = 0.1585, where tmax = 20. 
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^ Sim 






0.8 


0.472(3) 


0.502(3) 


0.030(1) 


1.0 


0.485(3) 


0.512(4) 


0.026(2) 


1.2 


0.494(3) 


0.517(4) 


0.023(2) 


1.7 


0.506(3) 


0.525(4) 


0.018(2) 


2.0 


0.510(3) 


0.527(4) 


0.016(1) 


2.5 


0.514(4) 


0.528(4) 


0.014(1) 


3.0 


0.517(4) 


0.529(4) 


0.012(1) 


3.5 


0.519(1) 


0.530(1) 


0.010(1) 


4.0 


0.520(4) 


0.530(4) 


0.009(1) 


7.0 


0.522(4) 


0.529(5) 


0.007(1) 


10.0 


0.522(4) 


0.528(5) 


0.005(1) 


oo 


0.529(6) 





Table 4: Ground state energies extracted from n^^p = 1 fits to Cn for = 
0.1585. Also shown are the hyperfine splittings obtained from performing a single 
exponential fit to the ratio of the ^Si and Cu correlators. 
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aMo 






a{2S-lS) 


1.0 


0.485(3) 


0.83(2) 


0.34(2) 


2.0 


0.510(3) 


0.81(2) 


0.30(2) 


4.0 


0.520(4) 


0.77(1) 


0.25(1) 



Table 5: Ground and first excited state ^Sq energies obtained from Uexp — 2 fits 
to Cii and C21 for ki — 0.1585. 
and first excited ^5*0 state. 



Also shown is the splitting between the ground 
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Figure 5: The effective masses of the ^5*0 meson correlators for aMo 
Hi = 0.1585. 
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Figure 6: Single exponential fits to the ^5*0 Cu meson correlator (circles) and 
simultaneously to the Cu and C21 correlators (crosses) for aMo = 00 and ki = 
0.1585. tmax = 11 for the fits to Cu, and 15 for the fits to Cu and C21. The 
corresponding values of Q are shown as a solid line for the fits to Cu and as a 
dotted line for the fits to Cu and C21. 



Cu and simultaneously to Cu and C21, detailed in figure ||. A Uexp = 2 correlated 
fit to the vector of smeared correlators with Q > 0.1 was not found. The ground 
state and first excited state smearing functions are not sufficiently different with 
these statistics at early timeslices to resolve the first excited state. Consistency is 
found between the vector fits and the single exponential fits to Cu and we choose 
the latter with the time range 3—11 as the best fit. 

The static case is the infinite quark mass limit of NRQCD. Thus, in the 
smooth transition from finite Mq to the static limit the signal /noise problems, 
associated with qMq = 00, will also present a problem at finite Mq when the 
heavy quark mass is large enough. However, as seen in figure |^ we found that 
even the highest value of Mq used is not hampered by noise and the signal for 
the ground state remained out to timeslices of 18-20. The static limit appears to 
be a special point with poor signal /noise properties and our conclusions about 
the static result will actually rely on the infinite mass extrapolation of the more 
accurate NRQCD results. 



5.2 Variation of Egim with heavy quark mass 

The results for Esim for both ^Sq and ^Si at finite Mq and in the static limit are 
plotted against I/Mq in figure |^. As Mq increases and the hyperfine interaction 
becomes smaller the S-states become degenerate and tend towards the static 
result where flavour and spin symmetry are restored. This agreement with the 
static result coupled with the multi-state, multi-smearing analysis at finite Mq 
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Figure 7: The effective masses of the ^5*0 meson correlators for qMq = 10.0 and 
Ki = 0.1585. 



suggests that the excited state contribution to Egim for aMo = oo is small. In 
addition, the deviation from the static limit around the B meson (oMq ~ 2.0 
see Section p.2.1| ) is small and of order a few percent as expected by naive power 
counting arguments. 

In the same way as for the NRQCD action, Egim and other meson quantities 
can be parameterised in the heavy quark limit in terms of an expansion in 1/M. 
The coefficients of the expansion are nonperturbative quantities which can be 
extracted from the simulation. From first order perturbation theory in [H^ + 6H) 
about the static limit we obtain 



E ■ = + 



'{P\Qi-D^-CBcr-B)Q\PY 



1 



2Mn 



(12) 



where |P)°° represents the meson in the limit of infinite heavy quark mass and the 
0{1/Mq) coefficient is merely the expectation value of the kinetic and hyperfine 
terms appearing in the NRQCD action. Note that the possibility of using a 



different choice of expansion parameter in equation ^ introduces an 0{1/Mq) 
uncertainty in the first order coefficient. -E^^ is related to the binding energy of 
the heavy-light meson. A, in the static limit: 



A = lim (Mftz-mg) 

rpco pco 



(13) 
(14) 



where uiq is the pole mass and is the energy of a heavy quark in the static 
theory. Both A and °°{P\ Q{—D^ — cbO' ■ B)Q |P)°° appear in expressions for 
0(1/M) corrections to HQET predictions, but cannot be calculated using HQET 
alone. 
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Figure 8: aEsim vs l/(aAfo) for ^5*0 (squares) and ^5*1 (diamonds) for ki = 0.1585. 
The static value is shown as a circle and the solid lines show quadratic fits to 
the heaviest six data points. The data points shown as crosses correspond to 
gMq = 3.5 and 4.0 and are not included in the fits. 



To extract these quantities we performed correlated fits to the data as a 
function of 1/Mq using the functional form 

= + ^ + § + ^^^^ 

The fitting procedure was similar to that for the propagator fits; beginning with 
a constant fit function the endpoint for the fit is fixed to the heaviest data point, 
aMo = 10.0, and the initial point is varied over all finite values of Mq. The 
procedure is repeated for a linear, quadratic and a cubic fit. The results for the 
^5*0 meson are shown in table ^, where, except for the fits to a constant, only 
'good' fits to the data are presented (defined as Q > 0.1). 

The variation of the data with I/Mq is consistent with a constant only for the 
heaviest two data points, while a linear term is required in the region of the Bg 
meson and a quadratic term for l/aMo^O.8. The lightest two points require a fit 
function of higher order than a cubic function to be included in a fit. Since we 
have truncated the NRQCD action at 0{1/Mq) the coefficients of the quadratic 
term are not correct and if the results are to connect with simulations around 
the D meson the 0(1/ Mq) terms in the action are needed. 

In general, in order to obtain stable values for Ci the fit must include Cj+i. 
The results for Cq from both the linear and quadratic fits are consistent, and also 
in agreement with the static result. There are not sufficient data points in the 
quadratic region to determine whether the value for Ci is stable around ~ —0.03. 
Hence, only a rough determination of the coefficients is possible and the results 
are detailed in table |^ for both and ^5*0. Quadratic fits to both -E^^ EYi^ 
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order 


fit range 


n.d.o.f 


Q 


aCo 




a^C2 





1-2 


1 


0.81 


0.522(3) 








1-3 


2 


0.04 


0.515(3) 






1 


1-3 


1 


0.55 


0.527(4) 


-0.029(8) 






1-4 


2 


0.61 


0.529(4) 


-0.036(5) 






1-5 


3 


0.29 


0.533(3) 


-0.044(3) 






1-6 


4 


0.44 


0.532(3) 


-0.043(2) 




2 


1-6 


3 


0.42 


0.528(4) 


-0.031(9) 


-0.018(7) 




1-7 


4 


0.32 


0.524(4) 


-0.018(6) 


-0.015(4) 



Table 6: Correlated fits to E^^^ as a function of I/Mq for ki = 0.1585. The 
fit range 1—3 denotes a fit to E^^^ at aMo = 10.0,7.0 and 3.0. The points at 
qMq = 4.0 and 3.5 are not included in the fits for reasons explained later. 





aCo 




Sim 


0.528(5) 


-0.03(1) 


E^. 

Sim 


0.528(5) 


-0.01(1) 


E ■ 


0.530(5) 


-0.02(1) 



Table 7: The coefficients extracted from fits to -Esim as a function of 1/Mq for 
Ki = 0.1585. The errors include the variation in the coefficients obtained using 
different orders in the fit function. 



are shown in figure ||. The spin-average of the energies, Esim = |(-E'iim + '^^Yim)y 
removes the spin dependence and the slope of Esim depends only on the kinetic 
energy of the heavy quark. The results for Egim are also given in table |^. 

It is useful to compare with previous results for the slope and intercept of 
Egim- However, most other calculations have been performed in the quenched 
approximation; it is more appropriate to compare our results at = 6.0 (us- 

ing the tree- level Clover action for the light quarks as detailed in reference 0]) 
with those of other groups. In particular, Hashimoto H] using NRQCD for the 



heavy quark and the Wilson action for the light quarks and Crisafulli et al |]T0 
using the static approximation and the tree-level Clover action have results at 
this (3 value. Note that both groups use naive operators i.e. do not use tadpole 
improvement. In addition, taking into account the different systematic errors 
inherent in each of the simulations only a fairly rough comparison can be made. 
Considering the intercept initially, 

aE^^ = 0.60^ Hashimoto (16) 

= 0.61(1) APE Collaboration 00 (17) 
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where ki = Kc- This is in good agreement with our results at the same (3 value of 
aE"^^ = 0.461(15) taking into account the difference of Iumq = —0.13 (where Uq = 
0.878 measured from the plaquette) since we use tadpole improved operators. 

In order to calculate A from this quantity we use the tadpole-improved per- 
turbative value of = 0.29(8) at /3"/=° = 6.0 0, g. E^ is an 0(a,) quantity 
and the corresponding error from omitting 2-loop corrections may be large, dom- 
inating the error in A. Here the error in E^ is estimated by assuming the 2-loop 
corrections are roughly the square of the first order terms. CrisafuUi et al sug- 
gest that as a linearly divergent quantity, nonperturbative effects may make a 
significant contribution to E^ and using a perturbative value will give rise to 
an ambiguity in the corresponding value for A. Using a nonperturbative renor- 
malisation procedure which aims to remove all mixing with lower dimensional 
operators they obtain the naive (non-tadpole improved) E^ = 0.521(6)(10) [p!0[] . 
However, the disagreement with the corresponding naive perturbative result, 
E^ = 0.29(8) — In Mo = 0.42(8), is only at the la level. In addition, the nonper- 
turbative result appears to be gauge dependent |T2[ . 



From equation we obtain, 

aA = 0.17(10) K = Kc = 6.0. (18) 

for our quenched results. The error is large enough to provide consistency with 
aA = 0.09(1) from CrisafuUi et al. 
For this study we find, 

aA = 0.17(14) K = 0.1585, 

= 0.13(14) K = Kc. (19) 

using E^ = 0.36(13) from perturbation theory. The results for the meson binding 
energy at finite values of Mq, Amq = Esim—Eo, are given in table || and will be 
referred to below. It is clear that the uncertainty in the perturbative value for Eq 
must be reduced before a meaningful determination of the meson binding energy 
is possible with this approach. 

Consider the slope of Egim with Mq: in reference we found results for Ci 
similar to those presented in table |^, 

a^Ci = a^ < -&/2 ^t-dpole-irnproved^ -0.02(1) /5"^=° = 6.0 (20) 

independent of ki. The slope is small and negative. The naive, bare kinetic 
energy operator, however, is positive definite; in mean-field theory the two matrix 
elements are related by 

^ n2 In tadpole— improved ^ f^2 /q ^naive o /-i „ \ /'oi ^ 



^Th.e error in this quantity is not given. 

^CrisafuUi et al use the APE results generated on the same configurations for this quantity. 
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aMo 


3.0 


3.5 


4.0 


7.0 


10.0 


oo 


Ki = 0.1585 


0.168 


0.171 


0.174 


0.172 


0.171 


0.172 


^/ — 


0.129 


0.133 


0.135 


0.132 


0.134 


0.131 



Table 8: oAmq for the heavier quark mass values. The corresponding errors are 



of the same magnitude as those in equation IS. 



Thus, we obtain an estimate for < —D'^/2 >harT^ +0.35. Hashimoto extracted 
the slope of Esim (the hyperfine term is omitted from the NRQCD action), at 
K = He, and found < —D'^/2 >harT— +0.33 in good agreement with our results. 
Furthermore, Crisafulli et al, using their nonperturbative procedure calculate the 
matrix element explicitly and find < —D'^/2 >1a7e — +0.36(7) for ki around 
Kg. The overall agreement between different groups is encouraging. However, the 
dominance of the tadpole contributions to < —D'^/2 >harT suggest it is better to 
quote the tadpole-improved matrix elements. 

The physical kinetic energy of the heavy quark within the B meson, 
< —D^ > renorm, is a positivc definite quantity and expected to be very small, 
of size 0{a'^AQ(jj-)) ~ 0.04. The renormalised matrix element corresponds to the 
slope of Amo with respect to Mq. From table |Tl|, Eq has a negative slope with 
Mq and is of a magnitude that indicates a positive result for < —D^ >renorm that 
is not inconsistent with 0(Aq^^). However, due to the large perturbative error 
our results for A^/y, shown in table ||, are consistent with a zero slope. 

Also shown in figure § are the results for E^-^ and E^-^ at aMo = 3.5 and 4.0. 
The stabilising parameter, n = 1, was used for both these values of Mq. While this 
value satisfies ti^S/Mq the results for Egim do not smoothly interpolate between 
the heavier data points with n = 1 and the lighter points which use n = 2. In 
particular, it was not possible to perform a correlated fit (even including higher 
powers in 1/Mq in the fit function) that includes oMq = 3.5 and 4.0 and the data 
at aMo = 3.0, 2.5, ... for which n = 2. The problem is remedied by increasing 
n for these values of aMo to 2 |[T^; the shift in Egim is much less than la in the 



statistical errors and it is only the high correlations between results at different 
Mq and our statistical accuracy that allows us to see this problem. Using n = 1 
is sufficient for the heaviest data points at qMq = 7.0 and 10.0 and a correlated 
fit to the data can be performed omitting aMo = 3.5 and 4.0. These data points 
are omitted from all further fits to quantities as a function of 1/M. 

5.2.1 Meson Mass 

The heavy- light meson mass, M2, can be calculated both perturbatively and 
nonperturbatively. It is important to find consistency between methods and 
investigate their range of validity, not least to confirm that Lorentz invariance 
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can be restored at this order through a mass shift to the heavy quark. We 
calculated M2 in three ways: 



1. Directly, from the dispersion relation. M2 is the mass which appears in the 
kinetic term of a non-relativistic meson, 

^(p1 = Mi + i^ + ..., (22) 

where Mi = Esim- For comparison with later results we define the heavy 
quark mass shift as A/j/ = M2—M1. We extract = E{j)) — Mi by 
fitting the ratio of finite momentum to zero momentum ^6*0 correlators 
for Cii and C21 simultaneously to multiple exponentials for |p| = 1 and 
\/2 (for a fit involving more than one exponential this is not quite the right 
ansatz; however, it should allow most of our excited state contamination 
to be removed from AE). The next order term in the dispersion relation, 
|p|^/8M|, is not correctly reproduced by our 0{1/Mq) action and thus 
M3 7^ M2. However, for lighter meson masses this term may be significant 
and we solve for M2 using 

assuming M2 — M3 is small. Errors on the meson mass are obtained simply 
by solving for M2 using AE ± Icr. 

As an example the results for ^5*0 at oMq = 1.0 and = 0.1585 for = 1 
are shown in figure |^. The effective mass illustrates a difficulty in ex- 
tracting M2 from the dispersion relation; the smearing functions which 
are adequate for zero momentum have a reduced overlap with the ground 
state at finite momentum |l^. The effective mass does not plateau for 
CiKIpI = = 0) until t ^ 12 (compared to t ^ 5 at zero- 

momentum), when the statistical errors are much larger than at earlier 
timeslices. However, the onset of the plateau is clear from a comparison 
of the effective mass using the excited state smeared source, which has a 
very different overlap with the ground and excited states. The variation of 
AE with tmin is stable from around tmin = 13 for nexp = 1 and consistency 
is found with the results obtained with additional exponentials. A fitting 
range of 5 — 20 with ng^p = 2 is chosen for oMq = 0.8 to 2.0 and 13 — 20 with 
^exp = 1 for aMo = 2.5 to 4.0; the corresponding values of AE are given 
in table |^. Finite momentum correlators were not generated for aMo = 7.0 
and 10.0. 

As the heavy quark mass increases, AE decreases, and correspondingly the 
determination of M2 is statistically more uncertain. Table |10| details the 
results for M2 from both = 1 and ^/2. Agreement between the meson 
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Figure 9: The effective masses of Cii{\p\ = 1)/Cii{\p\ = 0) (squares) and 
C2i{\p\ = 1)/C2i{\p\ = 0) (diamonds) for ^So at uMq = 1.0 and ki = 0.1585. The 
fit parameter AE{p) extracted from 1 (crosses), 2 (diamonds) and 3 (squares) 
simultaneous exponential fits and the corresponding values of Q are also shown. 
tmax is fixed to 20. 
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aMo 




fit range 


Q 


aAE{p) 


aAE'{p) 


0.8 


2 


5-20 


0.6 


0.060(2) 


0.37(7) 


1.0 


2 


5-20 


0.5 


0.051(2) 


0.33(6) 


1.2 


2 


5-20 


0.4 


0.045(2) 


0.30(5) 


1.7 


2 


5-20 


0.2 


0.034(2) 


0.26(5) 


2.0 


2 


5-20 


0.2 


0.029(2) 


0.25(5) 


2.5 


1 


13-20 


0.2 


0.025(2) 




3.0 


1 


13-20 


0.2 


0.021(2) 




3.5 


1 


13-20 


0.1 


0.019(2) 




4.0 


1 


13-20 


0.1 


0.016(2) 





Table 9: Fit parameters and values of AE{p) extracted from performing a si- 
multaneous correlated fit to the ratios Ci;(|p| = 1)/Ci;(|p| = 0) and 6*2; (|p| = 
1)/C2ii\p\ = 0) for Ki = 0.1585. 



mass from both momenta is found to be less than one standard deviation 



for aMo > 1-2, and confirms the dispersion relation in equation |2^. For 
the lighter meson masses the disagreement is less than 2a but indicates the 
increase in significance of the \p\'^/8M^ contribution and M2 7^ M3. Hence, 
we take M2 obtained from the lower value of momentum as more reliable. 
However, for meson masses around the B meson and less the statistical 
uncertainty in M2 is becoming prohibitively high. 

2. Perturbatively. Since the inequality of Mi and M2 is due to the omission of 
the constant term in the NRQCD action, Lorentz invariance can be restored 
by a shift to the heavy quark mass, 

M2 - Ml = Z^Mo - Eo. (24) 

The mass renormalisation, Zm, energy of a heavy quark in NRQCD, Eq, 
and the corresponding mass shifts, Ap^rt = M2 — Mi, as calculated in 



perturbation theory are shown in table |Tl| [^. The characteristic scales, 
and q^^, of Zm and Eq respectively are also given and these indicate the 
reliability of the perturbative result; if the characteristic scale is very soft 
perturbation theory is not a reliable method of calculating the quantity. For 
aMo around 2.0 and qMq = 7.0 and 10.0 the characteristic scales are hard 
and perturbation theory works well. However, as aMo 5, q^^ becomes 
softer. This defect is due to the method of calculating the characteristic 
scale and does not relate to any physical effect. We use the suggestion of 
Morningstar |F| to guess a characteristic scale q to smooth over the defect. 
The result A^^j is also shown in the table. We also used this method of 
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aMo 


aM2{\p\ = 1) 


aM2{\p\ = V2) 


0.8 


1.25(5) 


1.16(4) 


1.0 


1.49(6) 


1.40(4) 


1.2 


1.69(8) 


1.61(6) 


1.7 


2.25(14) 


2.20(10) 


2.0 


2.64(20) 


2.63(20) 


2.5 


3.07(27) 


3.19(29) 


3.0 


3.66(39) 


3.74(41) 


3.5 


4.05(48) 


4.27(53) 


4.0 


4.81(79) 


4.96(69) 



Table 10: M2 obtained from a\p\ = 1 and alp] = \/2 for ki = 0.1585. 



estimating the perturbative result for the lighter heavy quark masses where 
q* is soft because NRQCD breaks down. 



The mass shift is also shown in table |1^ along with an error estimating 
the higher loop corrections omitted in this calculation. The uncertainty 
associated with £"0 is computed as detailed in section |5l^ . For the quantity 
MoZm we estimate the two loop contributions to be roughly {Mq^Z^ — 1))^; 
this is a more conservative choice than the naive estimate of Mo{Zm — 1)^- 

From the dispersion relation of the heavy- heavy meson at the same Mq. 
As mentioned above M2 — Mi is due to a mass shift of the heavy quark 
and therefore the shift for a heavy-light meson should be half that for a 
heavy-heavy meson with the same bare heavy quark mass: 

Ahh = l[M2{hh)-Mi{hh)] (25) 

We used the NRQCD action including all relativistic corrections to 0(Mf*^) 
and the same method of extracting M2. The results are shown in table |12[ 
As mentioned in Section ^ discretisation errors in the heavy quark momen- 
tum become significant around aMo ~ 5 and thus the results for oMq = 7.0 
and 10.0 are not reliable. 



All three methods of computing the mass shift are compared in table |T2 
There is reasonable agreement for qMq < 4.0, while for the heavier quark masses 
agreement is not expected due to the large systematic errors associated with Ahh- 
The estimate of the perturbative result for aMo = 0.8 — 1.2 and 4.0 appears to 
work well. Furthermore, the agreement between obtaining the heavy-light meson 
mass directly and methods 2 and 3 which assume a mass shift can be applied to 
the heavy quark mass confirms that Lorentz invariance can be restored in this 
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aMo 


n 








aEo 




^ A"'' 

^'—^pert 


0.8 


4 


0.04 




0.03 






0.886 


1.0 


4 


0.19 


0.668 


0.14 


— 


_ 


1.063 


1.2 


3 


0.31 


1.755 


0.27 


1.481 


0.625 


1.254 


1.7 


2 


0.49 


1.253 


0.48 


0.402 


1.728 




2.0 


2 


0.54 


1.197 


0.55 


0.372 


2.022 




2.5 


2 


0.56 


1.146 


0.62 


0.359 


2.506 




3.0 


2 


0.52 


1.119 


0.68 


0.349 


3.008 




3.5 


1 


0.34 


1.155 


0.71 


0.348 


3.695 




4.0 


1 


0.23 




0.74 


0.346 




3.944 


7.0 


1 


25.10 


0.994 


0.79 


0.350 


6.608 




10.0 


1 


4.11 


0.978 


0.81 


0.351 


9.430 





Table 11: Perturbative results for M2 — Mi. A^g^^ is obtained using av{q*, q) 
where q = 0.6 and 0.8, and the resulting two values for the mass shift are averaged. 
For aMo = 0.8, 1.0 and 4.0, '-' indicates as is ill-defined at the corresponding 
value of q*z^ or q%^. 



way. Since using the heavy-heavy dispersion relation is the most accurate method 
we use Ahh for oMq < 4.0 and the perturbative results for the two heaviest quark 
masses. Using a^^ = 1.8—2.4 we find qMq = 2.4—1.7 corresponds to the bare b 
quark mass. 

5.3 Mass Splittings 
5.3.1 Hyperfine Splitting 

We extracted the hyperfine splitting by performing a single exponential fit to the 
jacknife ratio of the ^5*1 and ^Sq Cu correlators. The effective mass for qMq = 1.0 
and Ki = 0.1585 is shown in figure |1^ along with the hyperfine splitting extracted 
as train IS Varied. The plateau begins at t^j^^ll, and this is found to be the case 
for all oMq. The splittings extracted using a fitting range of 11 — 20 are shown 
in table |. 



Figure [Tl| presents the hyperfine splitting as a function of 1/Mps, where Mps 
is the ^Sq meson mass. Performing correlated fits to -E]^^ ~ ^sim using a fit 
function of the form in equation and the same fitting procedure we found 
linear dependence on 1/Mps around the B meson and down to oMq = 1.2. 
The extrapolation to the static limit is consistent, within 2a, with zero; the 
uncertainty in the meson mass has not been included in the fit. This is largest 
for the heaviest two points. If these are omitted the extrapolation to infinite 
heavy quark mass is closer to zero. The coefficients extracted were stable over 
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Figure 10: Shown are the effective mass of the ratio of ^5*1 to ^Sq Ch correlators 
for gMq = 1.0 and ki = 0.1585 and the mass sphtting extracted using Ug^p = 1 
as a function of tmin- The corresponding values of Q are also shown, tmax is fixed 
at 20. 
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aMo 


aAhh 




a^hi 


0.8 


0.89(2) 


0.89(34) 


0.78(5) 


1.0 


1.09(3) 


1.06(30) 


1.00(6) 


1.2 


1.27(2) 


1.25(31) 


1.20(8) 


1.7 


1.76(2) 


1.73(35) 


1.74(14) 


2.0 


2.07(2) 


2.02(30) 


2.13(20) 


2.5 


2.58(2) 


2.51(26) 


2.56(27) 


3.0 


3.10(3) 


3.01(25) 


3.14(39) 


3.5 


3.64(3) 


3.70(31) 


3.53(48) 


4.0 


4.12(4) 


3.94(20) 


4.29(79) 


7.0 


8.10(10) 


6.61(12) 




10.0 


13.80(25) 


9.43(22) 





Table 12: Heavy quark mass shifts determined from the heavy-hght dispersion 
relation, A/j^, perturbation theory, Apert and the heavy- heavy dispersion relation 



the fitting range and we find 

aCo = 0.016(7) a^Ci = 0.038(2). (26) 

where the coefficient of the slope corresponds to y < cbct ■ B >I^JI:{'°J:^^"^p^°'"^^ . 
This is consistent with OIo^Aq^^) ~ 0.04 expected from naive power counting. 
In addition, agreement is found with our results on quenched configurations, 
a^Ci = 0.043(6) Converting to physical units, 

<CBa-B >llt1^-'"'''''"'"'^ = 0.24 - 0.44 GeV' /?("^=2) = 5.6 (27) 

= 0.362 ± 50 ± 36 GeV^ /5W=o) = 6.0 (28) 

where for the quenched results the ffist error is statistical and the second is due 
to the uncertainty in a~^. This compares reasonably well with an explicit mea- 
surement of the naive matrix element obtained by Ewing et al [jl5l on quenched 
configurations using tree-level Clover fermions at ki = Kg and = 6.2: 

< a ■ B >= 0.465 1 25(stat) + gQ(syst) GeV^, where is the renormahsation 
factor required in the static theory. 

We find no dependence of the splitting on the light quark mass, and this is 
also found in experiment: B*-B = 46(1) MeV and B*-Bs = 47(3) MeV. We 
find a splitting of 25—40 MeV; the large uncertainty in the splitting reflects the 
particular sensitivity to the lattice spacing arising from the 1/Mps dependence. 
While the hyperfine splitting is uncertain it does appear to lie below the experi- 
mental value. This may be due to residual quenching effects, the correct number 
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Figure 11: The hyperfine splitting vs l/{aMps) for = 0.1585. The sohd hne 
indicates a hnear fit to the data shown as squares. The dashed hne indicates a 
hnear fit to the data (squares) omitting the heaviest two points. 



of quark loops being ~ 3 — 4. The splitting is proportional to the wavefunc- 
tion at the origin and is dominated by a harder physical scale than those which 
characterise the quantities used to calculate the lattice spacing. Hence, it is a 
quantity which is sensitive to quenching effects. 



5.3.2 P-S Splitting 

The first measurements of the P-states of the B mesons have recently been made 
at DELPHI and OPAL |T^, 0. This provides a further test of the lattice simula- 



tion of B mesons through comparison with experiment, and also a ideal quantity 
with which to extract a value for the lattice spacing from the B system itself; 
the P — S splitting is expected to be very weakly dependent on the heavy quark 
mass. So far the spin splittings of the P-states have not been clearly resolved 
and the experiments above find B**—B = 419(25) MeV, where a cross-section 
weighted mean is used for B**. A more detailed discussion of the comparison of 
our quenched results to experiment is given in 0. 

We obtain a signal for the ^Pi, ^Pq and ^Pi correlators; the ^P2 is too noisy to 
extract an estimate of the mass. However, the spin-splittings between the P-states 
are consistent with zero. This is shown in figure |l^ which displays the effective 
masses of the ^Pi Cu and Cu correlators and the ^Pi—^Pi and ^Pq—^Pi splittings 
for aMo = 1.0 and ki = 0.1585. The spin-orbit interaction of the light quark, 
which we expect to split the ji = 3/2 and ji = 1/2 states by 50 — 100 MeV |TB|, is 



below the level of our statistical errors. A multi-smearing analysis is not possible 
for the P-states, the first excited state smearing function used for the P states (see 
equation ^ had almost identical overlap with the ground and excited states as 
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Figure 12: The effective masses of the ^Pi Cu and Cu correlators, and the ratios 
of the ^Pi and ^Pi Cu correlators (squares), and ^Pq and ^Pi Cu correlators (di- 
amonds) . 



the ground state smearing function. 

As we are unable to resolve spin-splittings, we considered only the ^Pi state 
and performed Uexp = 1 fits to Cu and Cu separately. A simultaneous fit to 
both correlators was not possible since the Cu correlator becomes dominated by 
noise as Cu plateaus. The results, given in figure |1^, show a fairly significant 
discrepancy between the energies, stable in tmin, from the two fits. Unusually, the 
values from Cu lie above those for Cu, which does not approach a plateau from 
below. Since noise dominates the ^Pi Cu correlator around t ~ 10 it is possible 
that this correlator has not plateaued and the results from the Cu correlator, 
which can be checked over a larger range of tmin, are more reliable. 

In order to investigate the dependence of the P—S splitting on the heavy 
quark mass we calculated the splitting at aMg = 1.0, 1.7 and 4.0. From table |I3| 
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Figure 13: The ^Pi energy extracted from a Uexp = 1 fit to Cu (crosses) and 
Cii (circles) correlators separately and the corresponding values of Q (dashed and 
solid line for the fits to Cu and Cu respectively) for aMo = 1.0 and = 0.1585. 
tmax is fixed to 12 for Cu and 10 for Cu. 

the ^Pi — 5* sphttings, where S is the spin-averaged S-state, are 0{Aqcd) as ex- 
pected for the excitation of a light quark and only mildly dependent on the heavy 
quark mass. The results are plotted against 1/Mps in figure |T^. Considering the 
discrepancy of ~ 3a between the fits for Cu and Cu the dependence on Mq is not 
significant. We find no dependence of ^Pi — S on the light quark mass and using 
aMo = 1.7 we obtain a splitting of 400 — 700 MeV. The characteristic scale for the 
^Pi — S splitting is expected to be similar to that for light spectroscopy. However, 
using the experimental result to determine the lattice spacing, a"^ = 1.6(5) GeV, 
where the uncertainty due to fitting is included in the error. With such a large 
error this determination of the scale is not particularly useful and is consistent 
with all the other estimates of a~^. An improvement in the smearing functions 
and a resolution of the spin-splittings both in the simulation and in experiment 
is needed in order to provide a useful determination of the lattice spacing. 

5.3.3 Bs-Bd Splitting 

Another quantity of interest is the Bs—Bd splitting, which is expected to be very 
weakly dependent on the heavy quark mass at the 0{ms/M) level. In addition, 
as the splitting is the difference of the Bg and Bd binding energies it is expected 
to be 0(m,) ~ 100-300 MeV. Experimentally, Bg-Bd = 98(6) MeV. In order 
to extract the splitting we chirally extrapolated E^^^ using the energies from 
the data at ki = 0.1600 which are extracted using the n^xp = 1 fits to Cu 
correlators with the fitting range 3 — 20. Figure |15| presents our results for the 
splitting, which are consistent with no heavy quark mass dependence. Converting 
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aMo 


Cii 


Cn 


a(iPi - S) 


1.0 


0.74(1) 


0.81(2) 


0.28(6) 


1.7 


0.74(1) 


0.80(2) 


0.26(5) 


4.0 


0.73(1) 


0.80(2) 


0.25(5) 



Table 13: Energy of the ^Pi state obtained from n^j-p = 1 fits to Cu using the 
fitting range 5 — 15 and Cu using 3 — 10 for Ki = 0.1585. Also shown is the 
^Pi — S splitting, where the error bars take into account the uncertainty due to 
fitting. 
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Figure 14: The a(^Pi — S) splitting obtained using fits to the ^Pi Cu correlator 
vs l/{aMps) for ki = 0.1585. 
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Figure 15: The aMg—aMd splitting vs l/{aMps), where Ks = 0.1577 is obtained 
using the ratio M^/Mp. The burst indicates the static result which is not included 
in the fit. 



to physical units we obtain Bg—Bd = 90 — 130 MeV; consistent with experiment 
within the large systematic uncertainties arising from both and . 



6 Conclusions 

We presented a lattice study of the heavy-light spectrum in the region of the B 
meson, the first to partially include the effects of dynamical quarks. A multi- 
smearing, multi-state fitting analysis was performed and the ground S and P 
states as well as the excited S states were extracted. Figure |16| summarizes 
the results for the B spectrum. This is a high statistics calculation and the 
systematic errors arising from the Wilson light quarks, indicated in the figure, 
are the dominant uncertainty. Within these errors the results obtained are in 
agreement with experiment. We are in the process of repeating the calculation 
using tadpole-improved Clover light fermions. This will allow a comparison with 
our quenched results, in progress, which use the same light fermions. 

We performed a comprehensive calculation of the heavy-light meson mass, 
investigating three methods and their range of validity. The agreement found 
between these methods confirms that Lorentz invariance can be restored at this 
order in NRQCD by a constant shift to all masses. 

We calculated the spectrum over a wide range of heavy quark masses enabling 
a detailed investigation of heavy quark symmetry. For the binding energy and 
mass splittings the corrections to the heavy quark limit were found to be in agree- 
ment with naive expectations based on a hydrogen-like picture of the heavy-light 
meson. However, there is an indication from the mass dependence of the ground 
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Figure 16: The B spectrum. The filled circles denote our results, where = 
2.0 GeV has been used to convert to physical units and the error bars do not take 
into account uncertainties in a~^. The open circle indicates the estimate of the 
systematic uncertainty in the 2S — IS and ^Pi — S splitting. The dashed lines 
denote the upper and lower bounds on the experimental results. The B mass has 
been shifted upwards to match the physical value. 

S'-state energy that for quark masses around Mh the next order in the NRQCD 
expansion is required. Using our tadpole-improved approach the nonperturba- 
tive coefficients A and < —D^ > were extracted. We showed that previous lattice 
calculations have large tadpole contributions and when we correct for this using 
mean field theory we find agreement with our results. The major error remaining 
is a perturbative one. 
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